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(Dated: May 28, 2008) 

Abstract 

A precision measurement of the mass of the hd^Pi) state of charmonium has been made us- 
ing a sample of 24.5 million ip{2S) events produced in e+e^ annihilation at CESR. The reaction 
used was tp(2S) — > vr^/ic, 7r° 77, he "yr]c, and the reaction products were detected in the 
CLEO-c detector. Data have been analyzed both for the inclusive reaction and for the exclu- 
sive reactions in which r]c decays are reconstructed in fifteen hadronic decay channels. Consis- 
tent results are obtained in the two analyses. The averaged results of the present measurements 
are M{hc) = 3525.28 ± 0.19(stat) ± 0.12(syst) MeV, and B{ip{2S) ir^hc) x B{hc jr]c)= 
(4.19±0.32±0.45)xl0-'^. Using the ^Pj centroid mass, AM^filP) = {M{xcj)) - M(hc) = 
+0.02 lb 0.19 ±0.13 MeV. 
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The large body of experimental data for the spectroscopy of the charmonium (cc) states 
has provided detailed information about the QCD interactions between a quark and an 
antiquark. A convenient and transparent realization of the interaction is in terms of a 
potential which is generally assumed to consist of a Coulombic part attributed to a vector 
one gluon exchange, and a less well understood confinement part. In analogy with QED, 
the spin-dependence of the interaction is attributed to the Breit-Fermi reduction of the 
one-gluon vector exchange, which leads to spin-orbit (L ■ S), tensor (T) and spin-spin 
(Si ■ S2) potentials. The confinement part is generally assumed to be Lorentz scalar and 
no spin-spin dependence arises from it. The mass splitting of the triplet IP charmonium 
states into Xcoi^Po), Xci(^-Pi) and Xc2(^-P2) is determined by the (L ■ S) and (T) terms of 
the potential, and the (Si ■ S2) term determines the hyperfine or triplet-singlet splitting. 
If the qq hyperfine interaction receives no contribution from the confinement part, and is 
only due to the Coulombic term in the potential, it is a contact interaction in the lowest 
order, and it is identically zero for all L ^ 0, i.e., AMhf{lP) = M{^P) - M{}P) = 0. 
The triplet ^Pj states are well established, and the mass of their spin-weighted centroid is 
{M{^Pj)) = [M(xco) + 3M(xci) + 5M(xc2)]/9 = 3525.30 ± 0.04 MeV fl. The singlet state 
hci^Pi) was not identified until very recently [2I, [sj. Although the identification of the triplet 
centroid mass {Mi^Pj)) with the unperturbed triplet mass M{^P) has been questioned 
it is necessary to make a precision measurement of the mass of he irrespective of how M{^P) 
is determined. 

Two recent experiments have reported identification of he and measured its mass. The 
CLEO measurement [2| was made by means of the isospin-forbidden reaction 

^{2S) -> vfO/i,, 7r° ^ 77, he ^ 7^c (1) 

using 3 million ip{2S) produced in e^e~ annihilations. The he was identified as the enhance- 
ment in the mass spectrum of recoils against vr''. Two different kinds of analysis of the data 
were done. In the inclusive analyses he decays were identified by loose constraints on either 
the energy of the El photon from he decay, or the mass of rje- In the exclusive analysis 
no constraint was placed on -^(7). Instead, rje events were reconstructed in seven different 
hadronic decay channels of rje- The combined significance level of the he observation was >6 
cr, and the quoted mass was M(/;,c)=3524.4±0.6±0.4 MeV. 

The Fermilab E835 measurement made scans of antiproton energy for the reaction, 
pp ^ he 'yi]e, rje — > 77. The results from the year 1997 scan and the year 2000 scan were 
combined to obtain M(/;,c)=3525.8±0.2±0.2 MeV. The significance level of he observation 
was ~ 3(j. No evidence was found for he in the previously reported reaction pp —>■ he 

If it is assumed that M{^P) = {M{^Pj)), the above two measurements lead to 
AM/,/(lP)=+0.9±0.6±0.4 MeV (CLEO), and AM,,/(lP)=-0.5±0.2±0.2 MeV (FNAL). 
While both results are statistically consistent with the prediction of AM/j/(lP)=0, it is 
important to understand any deviation from it, and its origin. 

In this Letter we report a much improved measurement of the reaction in Eq. (1) using 
nearly an order of magnitude larger sample of N{ip{2S)) = 24.5 ± 0.5 million [6j obtained 
at the Cornell Electron Storage Ring with e~^e~ annihilations at a center of mass energy 
corresponding to the il^{2S) mass of 3686 MeV [l[. The CLEO-c detector was used for the 
detection of the reaction products. 

The CLEO-c detector [7| , which has a cylindrical geometry, consists of a Csl electromag- 
netic calorimeter, an inner vertex drift chamber, a central drift chamber, and a ring-imaging 
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FIG. 1: Spectra of recoil masses against 7r° in the inclusive analysis: (top) full spectrum, (bottom) 
background subtracted spectrum. 

Cherenkov (RICH) detector, inside a superconducting solenoid magnet with a 1.0 T mag- 
netic field. The detector has a total acceptance of 93% of An, photon energy resolutions 
of 2.2% at E^=l GeV, and 5% at 100 MeV, and charged particle momentum resolution of 
0.6% at 1 GeV. 

The event selection criteria common to both the inclusive and exclusive analyses are the 
following. The events were required to have at least three electromagnetic showers and two 
charged tracks meeting the standard CLEO quality and vertex criteria |8i]. The acceptance 
region was defined as |cos^| < 0.93, except that recoil 7r° candidates were reconstructed 
using photons only in the good barrel region, | cos^^| < 0.81. For showers it was required 
that E"^ (barrel) > 30 MeV, and i?^(endcaps) > 50 MeV, where the endcap region is defined 
as 0.85 < I cos 6^1 < 0.93. The events accepted for 77 decays of 7r° and rj were required to 
have M(77) within ±15 MeV of M(7r°)=135.0 MeV and M(r/)=547.5 MeV, respectively Q 
. It was further required that there be only one 7r° in the event with the recoil mass in the 
expected region of he mass, 3526 ± 30 MeV. These candidates were fit kinematically with 
M(77) constrained to the 7r° and 77 masses to improve energy resolution. To distinguish 
charged pions, kaons and protons a log-likelihood criterion including dE/dx and information 
from the RICH detector was used. 

In the inclusive analysis, in order to remove neutral pions from J /if) decays following 
ip{2S) Ti^n^J/ip and -k^ti^ J /tp, events were rejected with tt+tt" recoil mass in the range 
M{J/ip) = 3097±15 MeV and 7r°7r° recoil mass in the range M{J/i)) = 3097±40 MeV. Sim- 
ilarly, events with the invariant mass of all charged particles, M(all charged) > 3050 MeV, 
as well as events with recoil mass against 77 in the range M{J/ip) = 3097 ± 40 MeV, were 
rejected to remove decays through the xj states. 

For the inclusive analysis it is required that the energy of the El photon in he — »• jrjc be 
in the expected range -£^(7) = 503 ± 35 MeV. It is also required that there be only one such 
photon in the event. Further, this candidate photon was rejected if it made either a 7r° or 77 
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TABLE I: Results for the inclusive and exclusive analyses for the reaction ipi'^S) — > n^hc — > tt^^tJc- 
First errors are statistical, and the second errors are systematic. As described in the text, Bi x 82 
for the exclusive analysis is based on 165 it 19 counts. 





Inclusive 


Exclusive 


Counts 


1146 ± 118 


136 ± 14 


Significance 


IO.Oct 


13.2(7 


M{hc), MeV 


3525.35±0.23±0.15 


3525.21±0.27±0.14 


^1 X ^2 X 10^ 


4.22 ±0.44 ±0.52 


4.15 ±0.48 ±0.77 



with any other photon in the event. 

The mass spectra of tt^ recoils are shown in Fig. 1, with the full spectrum in the top 
panel, and the background subtracted spectrum in the bottom panel. When the requirement 
£'(7) = 503 ± 35 MeV is not imposed, a spectrum of the background is obtained with nearly 
twenty times larger yield and no apparent he enhancement, as is expected because of the 
small product branching fraction Bi{ip{2S) i^^h^ x -62 (/ic IVc) ~4xl0~^. To remove 
the small he contribution in the above background, events were removed if they had a 
photon with E{'~f) = 503 ± 50 MeV. In the fit of the he spectrum in Fig. 1 (top) this 
background shape was mapped to the full spectrum with just one normalization parameter. 
The peak shape used consists of a Breit-Wigner function with an assumed width of 0.9 MeV 
(same as r{xci)), convolved with the instrumental resolution function obtained by fitting 
the Monte Carlo (MC) simulation of the data. The xV^-O-f. of the fit is 54/52. In the MC 
simulations the angular distribution for the El photon was assumed to be (1 ± cos^ 6). The 
overall efficiency determined from the MC sample is e=ll.l%. The results of the fit, and 
BixB2 = N{he)/{t X N{iP{2S))) are listed in Table I. 

In the exclusive analysis no constraint on -£(7) was imposed. Instead, for the de- 
cays 4'{2S) n'^hc, he 'yrje, rje X, rjc candidates were reconstructed in fif- 
teen different decay modes, X, with multiplicities of 2 to 6. These modes were 
used because they had significant yields in the direct decays i'{2S) •yrjc. Sev- 
eral of them, marked with (*), were utilized for the first time. These decay channels 
are: pp, rjTi^Ti" [i] — > 77), r^vr+Tr^ (77 ti^tt^tt^), KsK^tt', K~^K~7r^, 7r+7r~7r"''7r~, 
K+i^"7r+7r", K^K^K^K~ , 71^71^71^71^71^11^ , 71^ 71^ 71^71^ , {*)T]K~^K^{r] 77), 

(*)pp7r°; (*)7r+7r"7r°7r°, (*)pp7r+7r", (*)7r+7r"7r+7r"7r°7r°. 

The decay chain in Eq. (1) as well as the above rje decays were identified from the recon- 
structed charged particles, and tt'^'s and rj^s. For rj decays to 7r+7r~7r°, it was required that 
the invariant mass be within 30 MeV of the nominal mass M(r/) =547.5 MeV |l|. For K'^ 
decaying into a 7r"^7r~ pair, it was required that the invariant mass of the pair be within 10 
MeV of the nominal mass M{Kg)=497.6 MeV [if, and information about vertex displace- 
ment was used to reject random tt+tt" combinations. The ip{2S) 71^71^ J/tp events were 
rejected with tt+tt" recoil mass in the range M{J/ip) = 3097 ± 15 MeV. 

The entire decay sequence was reconstructed for each r/e decay channel. A kinematically 
constrained fit was done for each event to take advantage of energy-momentum conservation, 
and it was required that the of the 4C fit be less than 15. The mass of the rje candidates 
was required to be within 30 MeV of the nominal mass of M('r]c)=2980 MeV |l|]. If multiple 
rjc candidates were found in an event, only the one with the smallest was retained. 
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FIG. 2: Summed distribution of recoil masses against vr*^ in the exclusive analysis with 15 decays 
channels of 7]c- See text for details. 

The 7r° recoil mass distribution for each decay channel was fitted separately using the 
instrumental resolution shape determined from MC simulation, convolved with a Breit- 
Wigner function of assumed width r(/ic)=0.9 MeV. The ARGUS shape ^ was used to 
parameterize the background. The fitted number of counts from individual decays range 
from 1 to 30. The summed distribution was fitted in the same way. The fit is shown in 
Fig. 2. 

The product branching ratio Bi{ip{'^S) ■n^hc) x B2{hc ■yrjc) is related to the observed 
counts in the different decay channels 77^ ^ X as the average 



N{X,h,)/e{X,K) \ 



Bi X B2 



^ N{X, direct)/e(X, direct) / ^ B{tlj{2S) 



(2) 



In order to minimize systematic errors in the evaluation of Eq. (2) it is desirable to 
construct rjc ^ X decays in the same manner for 7]^ from he and rjc from direct decay of 
ipi^S). We do so by placing a window of ±7 MeV around M{hc) in vr" recoil. The spectrum 
for the hadronic system mass for each individual decay channel was then reconstructed 
and fitted in the same manner for decays through he as for the direct decays. The fits 
were done using a Breit-Wigner function with r(?7c)=26.5 MeV convoluted with the 
experimental resolution function as determined by the MC simulation for that channel, 
and parametrized as a double Gaussian. The background in each case was parametrized 
using a polynomial. The number of counts in individual decays ranged from 37±11 (pp) 
to 1052±74 (7r+7r-7r+7r-7r°7r°), with a total SX(X, direct)=4043±127. The corresponding 
total SA^(X, he) was 165±19. This is larger than SA^(X, he) obtained by fitting the 7r° recoil 
spectrum for M{he) measurement, and has correspondingly larger efficiency. 

The efficiencies e{X,he) and e(X, direct) were determined from MC simulations sep- 
arately for each channel. As expected, it was found that the ratios of efficiencies, 
R{X) = e{X, direct)/e(X, he) were essentially independent of X, and had the average value 
(R) = 2.36 ± 0.17. This allows us to obtain from Eq. 2. 



Bi X B2 



EN{X, he 



B{ip{2S) ^ ^r^e) EiV(^, direct) 



X {R) = 0.096 ±0.013. 



(3) 



Using the summed counts above, and the recently measured CLEO value, B{ip{2S) 
7r/J = (4.32±0.67)xl0-3 we obtain B{^{2S) ir^he) x B{he jr]e 

(4.15±0.48(stat))xl0-^ 
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FIG. 3: Angular distributions of the photons from he — > 'jrjc. Circles and solid points denote results 
from inclusive and exclusive analyses, respectively. The curve shows A^(l + acos^9) distribution, 
corresponding to q=1.20, as explained in the text. 



TABLE II: Summary of estimated systematic errors and their sum in quadrature. N/A means not 
applicable. 



Systematic uncertainty in 


M(/le) 

Incl. 


, MeV B1XB2X 10^ 
Excl. Incl. Excl. 


Nms)) 


N/A 


N/A 


0.08 


N/A 




N/A 


N/A 


N/A 


0.66 


Background shape 


0.10 


0.01 


0.26 


0.15 


tt" energy calibration 


0.08 


0.08 


N/A 


N/A 


tt" signal shape 


0.03 


0.01 


0.14 


N/A 


he width 


0.03 


0.02 


0.27 


< 0.01 


efficiency 


N/A 


N/A 


0.20 


0.22 


Binning, fitting range 


0.03 


0.03 


0.08 


0.27 


MC input/output 


0.05 


0.11 


N/A 


N/A 


rjc decays 


N/A 


N/A 


0.18 


< 0.01 


r/c width 


N/A 


N/A 


0.16 


< 0.01 


rjc line shape 


N/A 


N/A 


N/A 


0.09 


Sum in quadrature 


±0.15 ±0.14 ±0.52 


±0.77 



The angular distributions of the El photons in both inclusive and exclusive analyses were 
obtained by fitting separately the he peak in the data for different angular ranges. The results 
are shown in Fig. 3. The distributions were fitted with the function A^(l + a cos^ 6). The fits 
give aind=0.S7±0.65 (xVdof=3.9/3) and aexd=l-89±0.94 (xVdof=1.8/3). In order to take 
the average of the results from inclusive and exclusive analyses, the exclusive events were 
removed from the inclusive sample. The average of the values from inclusive and exclusive 
analyses is Q!ai;erage=l-20±0.53, and the curve in Fig. 3 illustrates it. This is consistent with 
a=l expected for an El transition from hc{J^^ = V ) to i]c{J^'" = 0^+). 

Systematic errors in the two analyses due to various possible sources were estimated 
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by varying the parameters used. These include choice of background parameterization, 
r(^c) = 0.5 — 1.5 MeV, TT^ hne shape (varied by ±10%), bin size (varied between 0.5 
and 2 MeV), vr'' energy cahbration (varied energy of photon daughters by 0.2% to 1.0% 
depending on photon energy). In the branching ratio for '4'{2S) — 'jr]c (lO| the dominant 
systematic uncertainty is due to the hne shape of the rjc which propagates into our product 
branching fraction analysis. An additional 2% systematic uncertainty is included to account 
for the possibility that line shape for the El transition 'yrjc differs from that for the Ml 
transition in direct ipi^S) — >■ 7?7c in a way that does not cancel in Eq. (2). It was determined 
that the results were stable well within statistical errors for the variations of event selection 
criteria. 

The individual contributions to systematic errors, as well as their sum in quadrature, are 
listed in Table II. 

When the exclusive events are removed from the inclusive spectrum, and the data are 
refitted, we obtain M{hc) = 3525.35±0.27(stat). The average of this result for the (inclusive- 
exclusive) events and the result in Table 2 for the exclusive events gives our final result as 

M(/i,) = 3525.28 ± 0.19(stat) ± 0.12(syst) MeV, (4) 
Si(^(25)^7r%)xi32(/ie--7^c) 

= (4.19 ± 0.32 ± 0.45) x 10"^ (5) 

These results represent a large improvement over our earlier results. The significance of he 
identification is 10 a for the inclusive measurements, and 13 a for the exclusive measure- 
ments. The present results from the exclusive measurements are based on twice as many 
decay channels of t]c as before, and are in excellent agreement with the results from the 
inclusive measurements. 

The nearly one order of magnitude larger statistics available in the present measurements 
has enabled us to determine the systematic errors presented in Table II with much greater 
precision than in our earlier publication j^. This allows us to average the present results 
with the previous ones. The resulting average results are: 

M{K)ayg = 3525.20 ± 0.18 ± 0.12 MeV, (6) 
{Bi X i32)AVG = (4.16 ± 0.30 ± 0.37) x 10"^ (7) 

To put our results in perspective, we wish to make two further observations. 

It is expected that the El radiative transitions xa iJ/i^ and he — >■ 7?7c should be 
similar. Also, the total widths of T{xci) and T{hf^ should be similar. If we assume them to 
be identical, it follows that B2{hc —>■ 'jrjc) = B{xci iJ li') ~ 0.36 ± 0.02 fl[. Our product 
branching fraction then leads to Bi{il){2S) —>■ n^hc) ~ (1.13 ± 0.15) x 10~^. Incidentally, 
this is nearly equal to that for the only other isospin forbidden decay measured within the 
charmonium family, B{iIj{2S) —>■ 7r° J/-?/') = (1.26±0.13) x 10^^. A recent theoretical prediction 
(III gives the range B{ij{2S) 7T^he) = {0A-1.3)xl0-^. 

If the mass of the centroid of ^Pj states, {Mi^Pj)) is used as a measure of M('^P), the 
present measurement of M{hc) in Eq. 4 leads to 

AMhfilP) = (mCPj)) - MCPi) 

= +0.02 ± 0.19(stat) ± 0.13(syst) MeV. (8) 
The CLEO average mass in Eq. (6) leads to 

AMhfilP) = +0.08 ± 0.18(stat) ± 0.12(syst) MeV. (9) 
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These results are consistent with the lowest order expectation of IP hyperfine splitting 
being zero. We notice that the triplet mass used above was obtained as {M{^Pj)) = 
(M(3Po) + 3M(3Pi) + 5M(3p2))/9, which is the evaluation of M(^P) in the lowest order, 
when the spin-orbit splitting is perturbatively small. It has been pointed out 0] that with 
(M(3P2) - M(3Po)) ^ 140 MeV, the validity of the perturbative determination of M{^P) 
is questionable. Indeed, the perturbative prediction that M(^Pi) — M(^Pq) = |[M(^P2) — 
M(3Pi)] = 113.9 ± 0.3 MeV disagrees with the experimental resuh, 95.9 ± 0.4 MeV, by 
18 MeV. This necessarily implies that the true M(^P) is different from the centroid value 
(M(^Pj)). Since AM^fi^P) is expected to be small (~ few MeV), if not identically zero, it 
is important that higher order effects should be taken into account in deducing M(^P) from 
the known masses of ^Pj states [3], so that a true measure of AMhf{lP) can be obtained. 
Only then can the present measurement of M{hc) be used to distinguish between the differ- 
ent potential model calculations, whose predictions for AMhf{lP) vary over a large range 
because of the different assumptions they make about relativistic effects, the Lorentz nature 
of the confinement potential, and smearing of the spin-spin contact potential fl^]. Although 
the presently available lattice calculations do not have the required precision |13i], it may be 
expected that future unquenched lattice calculations will resolve these problems. 
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